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Algebraic random walks (ARW) and quantum mechanical random walks (QRW) 
are investigated and related. Based on minimal data provided by the underlying 
bialgebras of functions defined on e. g the real line R, the abelian finite group Zjv, 
and the canonical Heisenberg-Weyl algebra hw, and by introducing appropriate 
functionals on those algebras, examples of ARWs are constructed. These walks 
involve short and long range transition probabilities as in the case of R walk, 
bistochastic matrices as for the case of Zjv walk, or coherent state vectors as in 
the case of hw walk. The increase of classical entropy due to majorization order 
of those ARWs is shown, and further their corresponding evolution equations are 
obtained. Especially for the case of hw ARW, the diffusion limit of evolution 
equation leads to a quantum master equation for the density matrix of a boson 
system interacting with a bath of quantum oscillators prepared in squeezed vacuum 
state. A number of generalizations to other types of ARWs and some open problems 
are also stated. Next, QRWs are briefly presented together with some of their 
distinctive properties, such as their enhanced diffusion rates, and their behavior in 
respect to the relation of majorization to quantum entropy. Finally, the relation 
of ARWs to QRWs is investigated in terms of the theorem of unitary extension of 
completely positive trace preserving (CPTP) evolution maps by means of auxiliary 
vector spaces. It is applied to extend the CPTP step evolution map of a ARW 
for a quantum walker system into a unitary step evolution map for an associated 
QRW of a walker+quantum coin system. Examples and extensions are provided. 



* Quantum Probability and Infinite Dimensional Analysis: From Foundations to Appli- 
cations, QP-PQ Vol.18, eds. M. Schiirmann and U. Franz, (World Scientific, 2005), p. 
174-200. 
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1. Introduction 

Random walks formulated in an algebraic framework 1,2 ' 3,4 of finite groups, 
bialgebras and operator algebras as well as in the framework of Quantum 
Mechanics 5 — 13 , and references therein), arc investigated. Minimal data 
for such constructions consist of a bialgebra 14 and an integral (functional) 
defined on it, or alternatively of some Lie algebra, and two quantum systems 
modelling the walker and the coin system, together with a map modelling 
the coin tossing, that decides probabilistically the stepping of the walker. 

Examples of ARWs treated in the following subsections are walks on al- 
gebras of functions on R, on Zjv, and on the canonical algebra Hciscnbcrg- 
Weyl hw 15 (section 2). For those walks we show how to define the entropy 
functional of their respective integral and/or Markov transition operator, 
and how to deduce that these are entropy increasing random walks by using 
arguments based on the interrelations relations between majorization bis- 
tochastic matrices and entropy 16,17,18,19 . Moreover, a number theoretic de- 
composing of ARW on Z^r is analysed, that is called prime decomposition 20 , 
and refers to its factorization into products of similar smaller Zjv- walks. 
Mathematically this decomposition is based on the Chinese Remainder 
Theorem and the co-associativity property. Also for the ARWs on R, in 
addition to the usual case of short range walk with nearest neighbor (NN) 
transitions (Polya walk) 21 , we discuss in our algebraic framework the cases 
of: i) the NN centrally biased random walk (Gillis walk) 22 and the case of 
symmetric random walk with exponentially distributed steps (Linderberg- 
Shuler:LS walk) 23 . Finally, for the hw ARW, where its functional is con- 
structed by means of the eigenstates of the annihilation operator of that 
algebra i.e the family of coherent state vectors 24,25 , the continuous time, 
or diffusion like limit, is obtained 15 . This limit results into a trace pre- 
serving quantum master equation 26,27 for the density matrix of a quantum 
boson system, which physically is identified with the evolution equation of 
an open quantum boson system interacting coherently with a classical elec- 
tric field and incoherently (dissipative interaction) with a bath of quantum 
oscillators rigged initially into a squeezed vacuum state (squeezed white 
noise) 28,29,30 . 

Section 3, gives a concise prescription of the concept of QRW, using the 
example of QRW on integers as paradigm 13 . It briefly explains the notion 
of quantum coin system and the coin tossing map, and summarizes two 
emblematic properties of that walk, namely the quadratic enhancement of 
its diffusion rate due to quantum entanglement between the walker and 
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coin systems, and the entropy increase without majorization effect of its 
probability distributions (pd). This section ends with a group theoretical 
scheme of classification of various known QRWs. 

In section 4, a relation connecting ARW and QRW is put forward. The 
connection is grounded on the theorem due to Naimark that asserts the 
possibility of implementing in a unitary way a CPTP map operating on 
e.g the density operator of a quantum system 31 . This unitary extension is 
realized in the original vector space of the density operator augmented by 
an auxiliary vector space, the ancilla space in the terminology of Quantum 
Information theory 32 . Applied in the context of CPTP of a ARW, the 
ancilla space is identified with the quantum coin state space of a QRW. 
The Kraus generators determining the CPTP map of a ARW serve to built, 
albeit in non unique way, the unitary evolution operator of the associated 
QRW. This section concludes with the example of an explicit construction 
of a QRW associated to the hw ARW, of section 2. 

Finally, section 4, summarized some of the results and gives some 
prospected applications of the ARW-QRW concepts and formalism. 

2. Algebraic Random Walks 
2.1. The case of R 

Proposition 1. Let the bialgebra of real formal power series H = Fun(R) 
generated by the coordinate function X, and let the positive definite func- 
tional 4> '■ H — ► R, defined as <j) — ^2 iGZ Pi4> ai , with < Pi < 1, ^2 ieZ Pi = 
1, where = f{&i), the functional that evaluates any function 

/ e H, at the point on = ia, defined for some step a <G R+, for i e Z. The 
n-step convoluted functional becomes 

where p^ — D n ~ 1 p^ 1 \ with the stochastic column vectors p^ = 
(Pi k) )iez,k = 1,2,3,... and initially = p = (pi) ie z- Also D = 
(Dij) = (pi-j), Z a bistochastic infinite matrix (delta matrix), and 

4> = {4>ai)iez a column vector. Majorization ordering among pd's is valid at 
each step i.e -< , and consequently the ARW is entropy increasing, 

namely, S(^ n+1 ) > 5(0*"), S(T; +1 ) > 5(T«), where S(f n ) = S(T$) = 
S(p( n ^), with S(x) any Shur-convex function e.g the classical entropy i.e 

Proof: Operating with convoluted functionals on some function / e 
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^(33,34 35,36,37,38^9,40) yieldg 

ijez 

= E P^Pk-i^ = E P [ P D ^ k = 4> T Dp. (2) 
k,jez k,jez 

By induction we obtain the aimed relation <j>* n = <fi T p( n ) = <j) T D n ~ 1 p( 1 \ 
Two important properties of the delta matrix are: i) bistochasticity, i.e the 
column and row sums is one, which is expressed by means of the column 
vector of units e = (...,1,1,1,...), that is left and right eigenvector of 
D 1 i.e De = e, e T = e T D, and ii) the shift property i.e Dij = Pi-j = 
P»+i-(j+i) = This property in the case of ARW in Z governed by 

a pd with finite support, or in the case of a N C Z finite dimensional walk 
(see remark below), amounts to a bistochastic matrix D. 

The functional of the walk at each time step <j)* n = J2 ieZ p\ n>> (f>ai, is 
characterized by the pd = (p[ n ^) ie z , which in turn is determined 
by the bistochastic matrix D i.e p^ = Dp( n ~^ = D n ~ 1 p( 1 \ Let us 
assume that the pd's are of finite support (but see remark below), then 
by invocation of the theorem stating that two discrete pd's x = (x m ) m , 
y = (yk)k, that are connected by a bistochastic matrix D i.e x = Dy, are 
ordered by majorization 16 > 17 > 18 > 19 i. e x -< y, we conclude that the sequence 
of pd's {p( x \ p( 2 \ p( 3 \ . . .} , of site occupation probabilities resulting at 
each time step during the evolution of the walk is partially ordered by 
majorization i. e p^ -< p( n ~ x \ n = 1, 2, 3, . . . . 

Let us adopt now the definition of the entropy of a functional to be 
the entropy of the pd that determines that functional i.e we set S((f>* n ) = 
S(Tg) = S(p^), or more generally we do so for any convex function S : 
R — > R of the type of the so called Shur-convex functions e.g the classical 
Shannon entropy S(p) = J2iPi^°SPi> or * ne functions F(p) = J2iPi' f° r 
any constant k < 1, or F(p) = — YliPi 17 ■ 

By virtue of the theorem stating that x -< y implies F(x) < F(y) , 
where F(x) = J2if( x i)> f° r an y convex function / : R — > R, or other- 
wise said that the convex functions isotonic to majorization 17 , we conclude 
that for the pd resulting from the random walk the majorization ordering 
is valid at each step i.e. p^ -< p( n_1 ), n = 1,2,3, .. . , and this implies 
ordering for e.g their entropies i.e. S(p n ) > S(p n ~ 1 ), and similarly for their 
functionals and transition operators. As majorization order implies entropy 
increase, it is considered as a measure of disorder, and this allow us to con- 
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elude the ARW are getting more disordered in the course of time with re- 
spect to their site- visiting pd's, which are getting more entropic, approach- 
ing, if left uninterrupted, to the uniform distribution of maximal entropy. 
□ 

Remarks: 1) The assumption in the previous proof about the support of 
the pd's been finite is not actually necessary. In fact in the proof given by 
Hardy et. al 41 is stated that given p -< q for finite sequence of q, p pd's, use 
of Muirhead's algorithm leads to a bistochastic matrix A, such that q = Ap, 
and then the proposition: H{Ap) > H{p) for H information/entropy or 
more generally a Shur convex function, is applied. This proof is constructive 
and builds A in a number of steps not greater than the lengths of p, q, 
therefore for infinite pd's a theorem not using bistochastic operators for 
the characterization of majorization is needed. Such a theorem is provided 
in 42 . 

2) For the above proposition the corresponding Markov transition op- 
erator defined as — (<j> o A, is equal to = J2 ieZ Pie ai ~^ . For the 
simplest case of p± = (p, 1 — p), a± = ±a, and all other p's and a' s been 
zero, the continues limit lirrin^oa = (T^)™ = T£° has been obtained that 
leads to a diffusion equation 33 ' 3, 4,15 . 

3) The above general algebraic setting implies that the stepping prob- 
ability matrix D = (Dij) = (pi-j), i,j G Z, would be expanded in the 
enveloping algebra U(e(2)) of the Euclidean Lie algebra e(2) « iso(2), 
spanned by monomials of its generators {E + , E_, L}, that satisfy the defin- 
ing commutation relations 43 

[£+,£_] =0, [L,E±]=±E±. (3) 

An irreducible matrix representation of those generators in the Hilbcrt 
space 7i = ^(Z) spanned by the eigenvectors of the "distance operator" i, 
is useful in expressing the D matrix for various random walks, and look for 
solutions by means of e.g the Fourier method. This irrep in the canonical 
basis of TL 1 and using the same symbol for abstract generators and their 
matrices reads: 

L = rne m J m , E± = e m±1 e\ n . (4) 

Next we present three different random walks in R, that can be used 
as show cases of the scheme of ARW presented here. In these concrete 
examples the defining the walk transition probability matrix D 1 is written 
as an element of the U(e(2)) algebra. No attempt will be made to give an 
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algebraic solution for the problem of finding the n th— step site occupancy 
probability distribution, as this can be solved by other means. 

The examples include: i) the simplest case of symmetric nearest- 
neighbor (NN) random walk (Polya-walk 21 ); one of its deformations, ii) 
the NN centrally (site n = 0) biased random walk (Gillis-walk 22 ), which 
refers to a solvable case of a walk with no translational invariance, and 
stepping probabilities with a bias which has power law decay, or more 
specifically which decays in proportion from the origin of coordinates. The 
e— deformation parameter is chosen so that when s > 0, the walk is biased 
to enhance returns to the origin, while if e < 0, escape from the origin 
is enhanced, and iii) a symmetric random walk with non- nearest-neighbor 
transitions with transition step length decaying according to an exponential 
law (LS-walk 23 ). Explicitly we have: 

1) Polya-walk: symmetric nearest-neighbor(NN) random walk with 

Markov transition operator 

D P = \{E_+E + ) (5) 
with matrix elements the inter-site transition probabilities 

Dp(l,l') = ls iX _ 1 + ls iX+1 . (6) 

2) Gillis-walk: nearest-neighbor centrally (site n = 0) biased random 
walk with Markov transition operator 

D G = D G {e)= l -(E_+E + ) + (E_-E + )^P ( t, - 1< e < 1 (7) 

where Pq = e ej, P^ = 1 — P , the projection operators in the vector e , 
and its orthogonal subspace respectively in ^(Z), and with matrix elements 
the inter-site transition probabilities 

Dg{1^) = \{1 + y)5 iX _ 1 + \{1-j)5 iX+1 , ZVO (8) 

D G (Z,0) = ^,_i + ^,i- (9) 

3) LS-walk: symmetric random walk with exponentially distributed 
steps with Markov transition operator 

D LS = D LS (e) = e- k "{E k _+E k + ) (10) 

kez >0 
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with matrix elements the inter-site transition probabilities 



D LS (l,l') = ^- Y ^e 
= 0, 



l-l 
1 = 1 . 



> 



(11) 



Use of the eigenvector equations E±e = e, e' = e'E±, leads to the con- 
clusion that the transition matrices of the Polya, and LS random walks i.e 
the matrices Dp, Dp, and Dls respectively are bistochastic, while that of 
the Gillis walk D G , is column stochastic. Also from the definitions the fol- 
lowing two limits are deduced: lim e ^ D g (e) = lim^oo Dls(s) = Dp. 



Generalizations: The 2D generalized Gillis random walk with entangle- 
ment. This model describes a 2D NN random walk with a biased towards 
a point placed at (mi, m^) coordinates on the plane. The bias is an attrac- 
tive or repelling one depending on the sign of two parameters (£1,62) in 
reference to the motion along x, y axes respectively, and its strength de- 
cays following an inverse power law with characteristic exponents (ai,a 2 ) 
correspondingly. This is summarized by writing explicitly all parameters 
in the Markov transition operator D g d (ei, mi, cn;£2, m 2 ,a 2 ), and their 
domain of values {sj, rrij, aj} e { (—1, 1), Z, Z + }, j = 1, 2. The transition 
matrix along each axis is 

D G {e 3 , moi a ) = \{E_ + E + ) + (E. E+ )-^^—^P^ j = 1,2 

(12) 

where Pj = e^ej, Pj- = 1 — Pj, the projection operators in the vector ej, 
and its orthogonal subspace respectively in Z 2 (Z). The 2D transition op- 
erator consists of an entangled 32 convex combination of two factorizable 
ID transition operators with the parameters of position of bias site, char- 
acteristic decay exponent, and decay strength c.f {mj, aj, Ej, }, j = 1,2 
interchanged, it reads (0 < q < 1) 



D G = qD G (e 1 ,m 1 ,a 1 ) <g> D G {E2,m 2l a 2 ) 

+ (1 - q)D G (£2 7 m2,a 2 ) <g> D g (e\, mi, ai), (13) 

with matrix elements the inter-site transition probabilities 
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{l[ -mi)"> 



£2 



+(l-«) s 



1 + 



£2 



(/I - m 2 ) Q2 



(/a - m 2 ) a * 
1+ 61 



(/a -mi) Q i 



^ l , I ^ — 1 ^l"2 j £ 2 — 1 



^1 — 1^2,' 



(14) 



with (li,l 2 ) 7^ ( m i> m 2)i f° r the g component and (Z x , Z 2 ) 7^ (m2,mi), for 
the (1 — g) component of the convex combination. Also 

D cf{h,l 2 = m l,2'- h,l 2 = m 2,l) = + ^li,mi+l) 

•(^k.ma-l + 2^2,m 2 +l) + (1 - ^(g^i.ma-l + 2^i> TO 2 + l) 

■(2^2,mi-l + 2^2,mi+l)- (15) 

the matrix elements with (li,l 2 ) = (mi,m 2 ), for the q component and 
{li,l 2 ) = (m 2 ,mi), for the (1 — q) component of the convex combination. 
The role of bias and that of the entanglement of the two ID walks, can be 
investigated in an algebraic manner in terms of tensor product representa- 
tions of the iso(2) algebra, and it will be given elsewhere 44 . 



2.2. The case of Z N 

In this section we give a brief study of algebraic random walks on abelian 
groups Zn 3 , using their underlying bialgebra structure, and further inves- 
tigate possible forms of their decomposition into simpler and dimensionally 
lower ARWs, based on number theoretic properties of N . 

Proposition 2. Let the multiplicative abelian group Z^v = 
{e 1 g 1 g 2 ,...,g N - 1 } and the bialgebra H = C(Z N ) = Fun(Z N ) = C 
sparing 1 } . q 1 }, with dual algebra H* = CZ^ = Fun{Z* N ), with pairing 
given by evaluation. Let the positive definite functional (state) of a Z^v 
random walk = X)ie[o jv-i]P*5'i identified as weighted sum of elements 
of H, with < pi < 1, X^£[o N-i] Pi = 1- By means of the column vectors 
P = (j>i)ie[o,N-i],G = (3 l )i£[o,iv-i],we write it as <j> = G T P, where T de- 
notes transpose. Then the n-th step convolution becomes <j>* n = G T D n ^ 1 P, 
where D = (Dij) = (Pi-j\ mo d n), i,j € [0, N — 1] a circulant bistochastic 
matrix, to be called delta matrix. 



Proof: Straightforward (c.f. 3 ). 



February 1, 2008 10:36 Proceedings Trim Size: 9in x 6in 



RRRREVISED'Lcvico'Proc'ELLINAS 



9 



Remarks: 0) The delta matrix is more precisely a circulant bistochastic 
matrix 45 , that can be treated as an element of finite Hcisenbcrg groupiJjy 
20 , and this leads to an explicit solution for the dynamics of Zjv walk. 

1) Recall the following version of the Chinese Remainder Theorem 
(CRT) 46 : let N — N1N2, the decomposition of positive integer N, into 
product of coprimes Ni,N2, then working with the abelian additive groups 
of numbers Z Nl2 , mod7V 12 respectively, we can introduce the unique bi- 
jection 

5:Z N — >Z Nl xZ N2 , 5(x) = (x - pN u x - aN 2 ), (16) 

p, a G Z, that maps the numbers of Zn, into the ordered pair of its 
remainders after division by Ni,N 2 . Its dual map is the inverse p = 
<5 _1 = Zn 1 x Zm 2 — > Zn, which constructs the unique number x from 
its remainders (a, b), with respect to the divisors N\,N 2 , as p(a,b) = 
aNf Nl) + bNf N2) = x, where <p(c) = # |c > m e Z+- GCD(m, c) = 1| , 
the Euler function of given integer c, that equals the number of co-primes 
less or equal to c. 

2) The same factorization is valid for abelian multiplicative groups i.e 
Zn — Zjvi Zjv 2 , if N = N\N 2 , and N\,N 2 , are relative primes. 

3) Let for later use introduce now the map S — ► Vg, that uniquely 
determines from the CRT bijection 5, with S(i) — (ii,i 2 ), the isometric 
matrix V$ : C n — > C jvi <E> Cn 2 j written in the canonical basis as 20 

v s = J2 e m4 = J2 ^ ® e ^ e l 



and its inverse 



ii£Z Nl ,i?£Z N2 ii£Z Nl ,i2£Z N2 

(18) 

for which we have V 5 V} = V s V 5 -i = V 5oS -i = V idxid = l CjVl ® lc N2 , and 
VgVg = V s -iV s = V s -i oS = V id = l Cw • 

Notation: in the sequel and in order to distinguish the space dimen- 
sionality N, referring to certain c.g functional, operator, probability vector, 
co-multiplication etc, we will denote it by <j)[N],T[ N ], and P[n] respectively. 
We can now state two necessary and sufficient conditions, in order to ob- 
tain an isomorphic prime decomposition of a Zat ARW governed by a pd 
P[N], into a product of two others ARWs Zjv 17 Zjy 2 , with respective pd's 
P[n] — P[Nx] ® P[n 2 ]- The first condition is number theoretic and is about 
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the compositeness of the dimension number N of the probability distri- 
bution P[jv] that generates the ARW on Zjv, while the second condition 
is about its factorization into a tensor product of two pd's of appropriate 
dimensions. 

Proposition 3. i) Let N — N X N 2 be the prime factorization of a positive 
integer N, then if in addition to the isomorphism of abelian groups of Zn ~ 
Zn ± ®Zn 2 , we consider a probability distribution (pd) P[n] factorizable into 
a product of two others pd's P[jvi] > P[n 2 ] > namely such that VsP[n] = P[Ni] ® 
P[n 2 ] i or p[ N ^ — p\^^Pi^ 2 \ then the functional of a Zn algebraic random 
walk <f> [N] = J2 t ez N P l f f] 9[ N] , with < pf ] < l,£ i6Zjr pf ] = !> factorizes 
for every step n, namely <f>*^ w ^j^] ® 0[jv 2 ]i an(1 similar factorization is 
valid for the transition operator i.e Tj-^ w (8) Tj^y • 

Let Ni, N2 and A r 3 j be the co-prime factors of some positive integer 
N, then the decomposition Z^ w <g> Z^ 2 ® Zjv 3 ~ Z^ l ^ 2 ® Zjv 3 ~ 
® Zn 2 n 3 , is co-assosiative 14 as indicated in the last two equations, 
and if a pd P[jv] is considered for which the finest decomposition VsP[n] = 
P[n-i] ®P[N 2 ] ®P[N 3 ] > m terms of three others pd's is true, then the associated 
Zjv ARW is also decomposed at any step n, namely for its functional and 
transition operator respectively, the following co-associative factorizations 
are valid 

<f>[N] ~ 0[JVi] ® ^[JV 2 ] ® <£[JV 3 ] ~ 0[WiiV 2 ] ® 0[JV 3 ] ~ ^[iVi] ® ^[JV 2 ]V 3 ] 5 ( 19 ) 

Proof: i) The A = NiN 2 case: By means of the factorization property 
of the pd viz. 

iGZjy iGZn ii€iZiy 1 ,i2€zZN 2 

E W^.®^^^]®^], (20) 

ii£Z]y 1 ,i2^Z^ 2 

we deduce that 

i,j€.Z N 

~ E p f(i-j) e ^ e ^ ' e ii e t = ^fyi] ® D [iV 2 ] , (21) 
and therefore VsD" N ^Vg w -D^j ® ^[jv 2 ]' ^ or the delta matrix. This yields 
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the factorization of the functional 

0[JV] = G[n] P [n] = Gj N] v}VsP[ N ] 

« (V s Gf N] )Hv s P [N] ) = Gf Ni] t ® Gf N2] lP [Nl] ® P [Ay = [JVl] ® 0^2,) 

and similarly for the transition operator (with t{x ® y) = y ® x ) 

T[N] = {id®r <g> id) o {4>[n\ ® id[N]) &[n] 

~ (^[JVi] ® ^[iVx] ® </>[at 2 ] ® id[ Nl ]) o (A [Nl] <g> A [iVl ]) « T [JVl ] ® T [A r 2] . (23) 

it) The AT = N 1 N 2 N ? , case: Applying the CRT to N, in two different 
ways yields for the abelian group Zn two factorizations i.e Zn ~ Zn x ® 
Z W2 ® Zat 3 rj Z NlN2 <g> Ztv 3 ~ Zj\Ti ® Z N2 n 3 - This can equivalently be 
expressed by means of the bijection <5, as 

(<Wiv 2 ] ® «d) o ^[JViJVa.JVs] = (id ® ^[jv 2 ,jv 3 ]) ° 6[n u n 2 n 3 ]- (24) 
Performing the factorization of the pd P[ N ] , as before twice we have 

(^5[Ni,N 2 ] 

which implies that 

(Vs [Nl , N2] ®l)V S[NiN2tN3] Dlf N] ((Vs [Ni ^ 

(26) 

This factorization of the delta matrix leads to an equivalent factorizations 
of the functional at each time step of the Zn walk, i.e 

4>*[N] « P&N,] ® 0[jv 3 ] « ® 0[aUs] « ® ^ 2 ] ® 0*^] ( 2 7) 

Similar decompositions can be obtained for Markov transition operators. 
□ 

Remark: Interpretation of this result implies that the original random 
walk Zn, is isomorphically decomposed into the product of : i) three sim- 
ilar walks Zjy — ZjVi ® Zn 2 ® Zjv 3 , ii) or into the product of two walks 
of dimensions N1N2 and ./V3 i.e Zn 1 n 2 ® Zat 3 and iii) or into the prod- 
uct of two walks of dimensions Ni and N 2 N 3 i.e Zn 1 ® Zn 2 n 3 - As an 
example we take the Z 6 random walk on the vertices of a canonical 
hexagon the statistics of which is determined by a 6-dimension pd vec- 
tor Prgi. If this pd has been chosen so that there are two other pd's one 
3-dimensional P31 that generates a Z3 ARW on the vertices of a canoni- 
cal triangle, and one 2-dimcnsional Ppi that generates a Z2 walk on two 
points, such that P[ 6 ] = P[3]®Pp], then we can decompose the ARW 
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on the canonical hexagon as a product of two walks one on the trian- 
gle times one on the two-point set. Similar interpretations can be given 
to the prime decomposition of an ARW on a canonical polygon. E. 
g the decomposition Z 60 = Z 3 (g) Z 4 (g) Z 5 = Z i2 (g) Z 5 Z 3 (g)Z 2 o 20 . 

Problem: Let N = NiN 2 be the prime factorization of a positive integer 
N, for which the isomorphism Zn ~ Z^ 1 <g> Zn 2 , is valid, let us consider an 
Zjv ARW generated as previously by a pd P\m- It this pd is not factorizable 
into a product of two others pd's as in the remark above, but instead there 
are two pairs of pd's P[Ni]i P[n{\> an< ^ ^[-^2]' ^""[Ay sucn that the original 
pd is a convex combination of them i. e (0 < q < 1), 

V S P [N] = qP [Nl] ® P [N2] + (1 - q)P( Ni] ® P^, . (28) 

In the terminology of quantum information 32 , we have here two 
ARWs ZjVi,Zjv 2 , that are classically correlated (cc) and form a proba- 
bilistic decomposition of the ARW in Zjv, and write symbolically Zjy — 
gZ^ (g) Ztv 2 + (1 — I^Nt ® ^w 2 - A number of interesting problems arise 
in this context: the total dynamics and reduced dynamics of the compo- 
nents of the Zat walk; the problem of construction of measures of cor- 
relations among the components of Z^ walk; the problem of informa- 
tion(majorization) dynamics and information exchange among the pd's 
components of the Z N walk; the problem of asymptotics of the Z N walk. 

2.3. The case of hw- Algebra 

A bialgebra 14 A — A([i, 77, A, e) over a field k is a vector space equipped 
with an algebra structure with homomorphic associative product map /1 : 
A x A — > A, and a homomorphic unit map 77 : k — > A, that are related 
by [i o (77 ® id) = id = fi o (id ® r/), together with a coalgebra structure 
with a homomorphic coassociative coproduct map A : A — > .4 <£> .4 and a 
homomorphic counit map e : A — > fc, that are related between them by 
(e ® id) o A = id = (id0e)oA. Both products satisfy the compatibility 
condition of bialgebra i.e (fj, <g> /u) o (id ® r ® id) o (A ® A) = Ao /_t, where 
r(a; (g) y) = y <E> x stands for the twist map. If rj or e is not defined in A we 
speak about non unital or non counital Hopf algebra. 

Suppose we have a functional <fi : A — > C, defined on .4, let us define 
the operator : .4 — > .4 as T$ = (0 ® id) o A, then e o T$ = <fi, namely the 
counit aids to pass from the operator to its associated functional. From this 
relation we can define the convolution product ip * 4>, between functionals 
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as follows 



e o T^T,), 



e o (^) ® id) o A o (<p ® id) o A 

(0 (g) V) ® irf ® e) o (irf ® A) o A 

{(f>®^})o/S. = (j)*lj), 



(29) 



and in general e o T£ = e o — <p* n . These last relations imply that 
the transition operators form a discrete semigroup with respect to their 
composition with identity element T e = id (due to the axioms of bialgebra) 
and generator T^, while the functionals form a dual semigroup with respect 
to the convolution with identity element e and generator <p, and that these 
two semigroups are homomorphic to each other. 

We recall now the Heisenberg- Weyl algebra hw and its structural maps: 
this is the algebra of the quantum mechanical oscillator and is generated 
by the creation, annihilation and the unit operator {a\a 1 1} respectively 
which satisfy the commutation relation (Lie bracket) [a, a'] = 1, while 1 
commutes with the other elements. This algebra possesses a non counital 
bialgebra structure 3 , chapt. 3), with comultiplication defined as 

A (n_1) a = n~5( a( g) ...<g>l + l(g>a<g>---<g>l + l<g>---<g>a), 
A^-^at = n -5( a t (g)---(gil + l(g)a(g)---(g)l + l(g)---(g)a t ), 



where as indicated above the A^™ -1 ) maps the creation/annihilation op- 
erators into the nth fold tensor product of the algebra and adds appro- 
priate factors (also c.f. 47 ). Let us also define the number operator 
N = a^a with the following commutation relations with the generators 
of hw: [N, a^] = , [N, a] = —a. The module which carries the unique 
irreducible and infinite dimensional representation of the oscillator algebra 
is the Hilbert-Fock space Hf which is generated by a lowest (or "vacuum" 
) state vector |0) € Hf and is given as Hf = {\n) = (a J |0) ,n € Z + }. 

The functionals we intend to use will be defined by means of the canon- 
ical coherent state vectors of the hw algebra so in the sequel we give a 
brief introduction to the concept of coherent state vectors (CSV) on Lie 
groups: consider a Lie group Q, with a unitary irreducible representation 
T{g), 9 € G, in a Hilbert space H. We select a reference vector |^o) € H, 
to be called the "vacuum" state vector, and let Go C G be its isotropy 
subgroup, i.c for h € Go, T(h) |*o) = e i{p{h) |*o>- The map from the fac- 
tor group M. = G/Go to the Hilbert space H, introduced in the form of 
an orbit of the vacuum state under a factor group element, defines a CSV 



Al 



1® 1 , 



(30) 



February 1, 2008 10:36 Proceedings Trim Size: 9in x 6in 



RRRREVISED'Lcvico'Proc'ELLINAS 



14 



\x) = T(Q/Qo) l^o) labelled by points x £ M. of the coherent state mani- 
fold. Coherent states form an (over)complete set of states, since by means 
of the Haar invariant measure of the group Q viz. dp.(x), x £ M., they 
provide a resolution of unity, 1 = J M dp(x) \x)(x\. As a consequence, any 
vector \if?) £ Tt is analyzed in the CS basis, |^) = J M dp(x)^(x) \x), with 
coefficients ty(x) = (a;|*). Wc should note here that the square integrabil- 
ity of the vectors "J/ will impose some limits on the growth parameters of 
the functions &(x) at the boundary of manifold M (cf. 24,25 and references 
therein). 

The hw-CS is defined by the relation 



oo 



\a) = e aa '~ &a |0) = Me aa " |0) = e^M 2 V 2= \n) . (31) 

It is an (over)complete set of normalized states with respect to the measure 
dfj,(a) = ±e-W 2 d 2 a for the non-normalized CS, and a £ M = HW/U{1) w 
C is the CS manifold. Since a\a) — a\a), M. is the flat canonical phase 
plane with the standard line element ds 2 — dada. Also the symplectic 
2-form cj = idahda is associated to the canonical Poisson bracket {f,g} = 

ndfdz - did*). 

The density operator (state) p which would be used to determine func- 
tionals of some operator bialgebras A, is defined generally as follows : Let 
a Hilbert vector space H that carries a unitary irreducible representation 
of A of finite or infinite dimension. The set of density operators 

S = {p£ End(H) : p > 0, p f = p, trp = 1} , (32) 

namely the set of non-negative, Hermitian, trace-one operators acting on 
"H form a convex subspace of End(W), which is the convex hull of the set 

S P = { P £S,p 2 = p} = H/U(l), (33) 

namely of the set of pure density operators (states), that are in one-to-one 
correspondence with the state vectors of H. 

The density operator to be used in the case of hw walk uses the pure density 
operators |±a)( ± a\ £ Sp and is a convex combination belonging to the 
convex hull of Sp i.e < p < 1, 

p = p\a)(a\ + (l-p)\-a)(-a\ . (34) 

Let </>(■) = Trp(-) =< p, ■ >, a functional defined on the enveloping 
Hciscnbcrg-Weyl algebra U{hw), where p = p\a){a\ + (1 — p) \—a){ — a\, 
i.c the p density operator is given as a convex sum of pure state density 
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operators. The action of the transition operator = ((f) ® id) o A on the 
generating monomials of U(hw) (where we ignore the numerical factors in 
the comultiplication of eq.(30)) reads, 

T ((o t ) m a n ) = [<t> <8> id) o A((a t ) m a") 

m\ / n 



[pa* 1 ^ + {l-p)(-a) l (-ay}{al) 



= EE 

= p(a t + a*)" 1 (a + a)" + (1 - p)(at - a*) m (a - a) 



1\m-i a n-j 



(35) 



For a general element f(a,a^) € U(hw) that is normally ordered, namely 
the annihilation operator a is placed to the right of the creation operator 
a\ denoted by f(a,a^) = J2 m n>o c mn(a^) m a n , the action of the linear 
operator becomes 

T (/(a, at)) = pf(a + a,a) + a*) + (1 - p)f(a - a, a* - a*) (36) 

By means of the CS eigenvector property and the normal ordering of the / 
element we also compute the value of functional viz. 



<t>{.f{a,a))) =pf(a,a*) + (1 -p)f(-a,-a*) . 



(37) 

Let us consider the displacement operator D a = e aa '- a ~ a = e A ; with 
A = aa) — oia which acts with the group adjoint action on any element / 
of the U(hw) algebra viz. 24 

AdD a {f) = Ade aaf - a * a (f) = e ««i(aat-a'a)(/) = e adA = Da f D ^ ( 38 ) 

where ad(X)f = [XJ] and ad(X)ad(X)f = [X, [X,f]] and similarly for 
higher powers, stands for the Lie algebra adjoint action that is defined in 
terms of the Lie commutator. Similarly the group adjoint action in terms of 
the displacement operator on the generators of U(hw) reads AdD± a (a) = 
a=F a and AdD± a (a') = a' =F a*. By means of these expressions we rewrite 
the action of the transition operator as 

T (/(a, a f )) = \pAdD- a + (1 - p)AdD a ]f(a, a*) 
= [pe- adA + (1 - p)e adA ]f(a, a*). 

Next we compute the limiting transition operator 

T t = Ts t = lim Tl l = 



(39) 



lim 

n— >oo 



1 + -ad A + -7 (acta 1 ) 2 + - 7 (acta) 2 - 
n n n 



t 



\j\ (ada ada^ + ada^ ada^) . (40) 
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In the last expression we have introduced the parameters of continuous time 
ieR and the drift and diffusion terms respectively c, 7 € C by means of 
the relations 

„ , 1. tc a 2 to , , „ , 

2 ^-2 )= n'T = ^ (41) 

and have performed the limits a — > 0, n — > 00, with t, c, 7 been fixed. 

We have also been used the limit lim Tl _ >0O (l + ff) = e z , to obtain the 

continuous time Markov transition operator T t = e tL , with its L generators 

as it is obviously identified in the equation below 

T t = e tL = cxpi [1 + adA + -f(ada^) 2 + j(ada) 2 - \ j\ (ada ado) + ado) ado))] . 

(42) 

By construction T t is the time evolution operator for any element / of 
U(hw) i.e f t = T t (f) and forms a continuous semigroup T t T t > = T t +t' under 
composition. This yields the diffusion equation obeyed by ft , which will be 
taken to be normally ordered hereafter. By time derivation of the equation 

Mf) =< P, ft >=< P, e tadC f >=< e- tadc 'p, f >=< p t ,f > , (43) 

we obtain the diffusion equation J^/t = £/t, as well as the dual one satisfied 
by the p density operator viz. ^p t — C) p t - 

Explicitly the quantum master evolution equation for the density matrix 
reads 

/!?(£) = [ca) — ca, p] + 7(a^ 2 p + pa) — 2a) pa)) + j(a 2 p + pa 2 — 2apa) 

- | 7 | ((2 AT + 1) p + p{2N + 1) - 2a) pa - 2apa)). (44) 

The obtained equation is similar to the quantum master equation that 
describes the trace preserving dynamics of the reduced density matrix oper- 
ator of single mode of the electromagnetic field interacting coherently with 
classical electric filed while it is immersed in a bath of quantum oscillators 28 . 
The decay of the field mode is influenced by the kind of initial condition 
the reservoir oscillators are put in. To analyze the physical content of that 
equation we rewrite it below by separating its right hand side into three 
lines i.e 

-^p(t) = [ca) - ca, p] 

— I7I {a) ap + pa) a — 2apa) + aa) p + paa) — 2a) pa) 

+ 7(a t2 /? + pa) - 2a) pa)) + j(a 2 p + pa 2 - 2apa). (45) 
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The first line gives the coherent interaction of the mode with the classical 
electric filed of intensity c, as described by the commutator of density op- 
erator with the Hamiltonian term. It is neglected for balanced walk. The 
second line is a typical part of a master equation describing mode decaying 
for reservoir oscillators in thermal equilibrium 48 . The last line is related 
to the case where the reservoir is prepared in a squeezed vacuum state 28 . 
Closing we should notice that the above quantum master equation can 
be transformed into a Fokker-Planck partial differential equation for some 
quasi-probability function e.g P, Q, or Wigner function associated with the 
density operator e.g 28 . 

3. Quantum Random Walks 
3.1. The case of Z 

In a quantum random walk there are two dynamically coupled systems: the 
walker systems described by a Hilbert space H w , and the coin system also 
described by a 2D Hilbert space H c ~ C 2 =span(|+ >, |— >). Let the uni- 
tary matrix U operating in H c , e.g Uh = ^ J ^j,or!7» = cf^Uh, the 

Hadamard (Fourier) transform and the j— rotation matrix respectively. 
If we denote by P± = |± >< =f|, the projection of an orthogonal partition 
of H c , and by S T , two step operators in the walker's space H w (explicit 
examples determined below), we introduce the unitary one-step evolution 
operator acting in the space H c <g> H w in the combined coin+walker system: 



In the equation above the upper/lower signs correspond to the choices 
Uh, U„ , respectively. If initially the two systems are decoupled, their den- 

4 

sity matrices are factorized i.e g c ®g w . Then if we assume that p c = |* >< 
^|, is a projective density matrix with |W >= a\+ > +b\— >, a normalized 
coin state vector, then the 1-time step of the QRW is considered as a com- 
pletely positive trace (CPTP) map e v , operating on the walkers' density 
operator, obtained by partially tracing out ("forgetting") the coin system 
i.e 




(46) 




(47) 



m— ± 
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Unitarity of V, implies that the two probabilities p± = TrP±U p c W P± = 
i ± 2Rca6 > 0, p+ + p- = 1, are determined by the coin system 
state vector variables and in turn they determine the Kraus generator 31 , 
(y/p+S + , y/p~ZS-), of CPTP evolution map. Four sources of choices are im- 
plicit in the above prescription of QRW: the choice of the initial coin+ walker 
state vectors, the choice of unitary U, the choice of definition of time step 
in terms of the tracing of the coin system (to be investigated in detail be- 
low), and finally the choice of step operators S T in the walker's space, that 
determines the king of the QRW under investigation. 

Next we choose to turn to a Hadamard random walk on integers with 
dynamical algebra the Euclidean algebra e(2) « iso(2), with step oper- 
ators ( v /p+S'+, ^JpZS-) = (y^P-E+j -\/(l — p)i?_), and third element the 
"distance" operator L. The latter is the interesting quantum observable 
the quantum moments of which are used to compare classical and quantum 
walks. We examine three possible tracing schemes: i) the classical scheme 
that promptly traces the coin system after each V action, leading to the 
CRW 

e v ( p W) = Tr c (V Pc ® p^V*) = pE +e ^E\ + (1 - p)E_ e ^E^ 

(48) 

which produces the diagonal sequence of density matrices 

{^),ev(pl 0) ),4(pl 0) ),4(pl 0) ),4(pl 0) ),---}, (49) 

with diagonal elements the probabilities of site occupancy, given by the 
(p, l—p) classical Pascal triangle, ii) the scheme that traces the coin system 
by increasing delays i.e after an increasing number of actions of V operator, 
leading to a walk designated by QRWi 

e VN (p^) = Tr c (V N Pc ® pEM") = £ fiW^t, (50 ) 

which produces the sequence of density matrices 

{pL 0) , e v (p^), e v2 (/i 0) ), e v >(pM), e v *{pM), ...} , (51) 

iii) and the scheme of delaying the trace of coin system by exactly one action 
of V operator, leading to a walk designated by QRW2 

e v -4pL n) )=Tr c (V 2 p c ®p[:- 1 ^) = B + g^Bl+B^^Bi, (52) 
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which produces the sequence of density matrices 

{ p<?> , e v , (pL 0) ) , 4, (p<?> ) ,4- ) , eU (Pl 0) ) , ■ ■ ■ } • (53) 
Next proposition deals with the time evolution of the pd (Pq N ^) m = 

(m\ | to) , made of the diagonal elements of the density matrices in the 
course of the walk QRWi and QRW 2 13 . 

Remarks: 1) From the above treatment can been shown that the QRWi 
reproduces along the diagonal elements of the sequence of evolving density 
matrices of the walker system the pd of the ID Hadamard random walk and 
its diffusion rate 13 2) by expressing the evolution unitary operator V, as 
the exponential of a hermitian operator H i.e V — e lH , which describes an 
interaction between coin and walker quantum systems, the various schemes 
of partial tracing can be physically implemented by choosing the length of 
interaction time e.g the QRW2 scheme requires an interaction time t = 2 
for the evolution operator V t = cxp [itH). AA 

Proposition 4- There exists bistochastic matrices A q = So o So + Si o Si , 
and Aq = Bo 0B0+B10B1, which determine the pd drawn from the diagonal 
elements of the evolving density matrices of the walks QRWi , and QRW2 
respectively, by means of the respective equations 

p(n+l) = Ac p(n) + ( E+ _ E_)MWPW , (54) 

P,(" +1 >=A,PW, (55) 

where 

M(") = V^T^ysW o S[ N) + S[ n) o S n) + (2p 1)S^ o S n) , (56) 
and 



A 9 = 



pa + \/p(l -p)b 



El + 



\/p(l - p)a - pb 



+ (1 - p)a - y/p(l - p)b + (1 -p)b+ y/p(l -p)a 1 



(57) 



Above the element by element or Hadamard product A = A o B, defined 
between matrices A, B of the same size by (^4 o B)ij = AijBij, has been 
used 17 . The study of the pd obtained from the previous proposition re- 
veals two novel aspects of the models QRWi and QRW2 . First the aspect 
of breaking the condition majorization-implies- entropy increase, and that 
of enhanced diffusion rates. Before closing this section we give a brief 
demonstration of the latter one (detail investigation together with relevant 
references can be found in 13 ). 
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Enhanced Diffusion Rates:Let us consider the mth order statistical mo- 
ment of the distance operator < L m >„= Tr{pffl L m ) , at the nth step. As- 
sume we have a symmetric walk with < L >„= 0, so that the standard devi- 
ation at nth step is a n = V< L 2 >„. For the CRW we have that = y/n, 
and from the pd's taken by the previous proposition we obtain the re- 
spective standard deviation for the walks QRWi and QRW2, expressed 
in terms of their classical counterparts, for the first five steps:{<7^ flH/l = 

}, and {*? RW > = a?,aQ RW > = ,a$ RW ° = y^g,**™' = 

j7]2o2,o-<i RW i=2c-C}. 

From that we deduce that there is a quadratic speed up of the spreading 
rate in the case of QRWi with respect to CRW, and that rate is even bigger 
for the case of QRW2. Finally the asymptotic (n 3> 1), growth values are 

turn out to be a^ RWl ~ ^Jn{2~V2)/2 of. 

4. Relation of Algebraic and Quantum Random Walks 

This section will put forward a relation between the two types of random 
walks under investigation so far i.e the ARQ and QRW. The relevant theory 
here is Naimark's extension theorem that allows to express in a non unique 
manner a positive trace preserving map, operating by means of its Kraus 
generators on a density matrix describing the state of some quantum system 
in a certain Hilbert space, by a unitary operator acting on a extension of 
the original space. 

Stated in the language of random walks the extension theorem assumes 
a ARW described by a CPTP map e(g w ) — J2m=±PmS m QwSl l , operat- 
ing on the density matrix of walker system with its Kraus generators 
(y/p+S + , y/pZS-), defined to act on Hilbert space H w . It is further as- 
sumed as usually for ARWs, that the step generators S± are related to and 
algebra of operators that needs to be specified. Then a unitary operator V 
is considered acting on H c Cg) H w i.e an extension of the original space by an 
extra or ancilla space H c C 2 =span(|+ >, |— >), which in the context 
of QRW stands for the coin system. Let a pure density matrix in the coin 
system p c — |W >< Then the extension theorem provides a unitary 
representation of the CPTP i.e 

£v{q w ) = ^ PmSmQwS^ = Tr c (V g c <g> g w V^). (58) 

The unitary operator provides the Kraus generators as y/p m S m — 
(m\ V |^) , up to a local unitary operator W, i.e the transformation V — > 
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W<g)l w V, provides the same generators. For the case of ARQs the unitary 
operator V is specifically expressed by means of the coin states projections 

P±, and the unitary matrix U(p±) = ( ^+ Vp^_ J of the coin space, as 



V = E m= ± PmU ® S m . In particular the step operators are unitary and 
inverse to each other i.e 5+ = (S-y . This unitary representation can be 
extended to products of CPTM maps by first defining the unitaries 

V® 2 = J2 P m U ® PnU ® S m S n . (59) 

m,n— ± 

Then we obtain 



Ey(ft»)= P™P n S mS n Qw(S m S n ^ =Tr c ®Tr c {V®Vg c ®Q c ®g w V^®V f ). 

m,n— ± 

(60) 

This unitary extension of e 2 , requires a double ancilla space or two coin 
quantum systems coupled, so the total space is H C ®H C ®H W . In the general 
case the fcth power of the CPTP map of an ARW can be implemented uni- 
tarily by extending the original walker space by fc anchillary coin systems, 
so the total space becomes Hf k £g> H w , and the total unitary operator is 

S rni •••S m , k (61) 

mi , . . .,mfc=i 

with unitarity condition V® k V® k t = lf k O l w . Then we obtain for the fcth 
step of the QRW as described by k successive actions of its CPTP map, a 
unitary realization which involves tensoring of quantum walker to k quan- 
tum coin systems, followed by a coupling of them by a unitary evolution 
operator on the space of coins+walker composite system, and finally a de- 
coupling of coins from the walker system, taken by partially tracing with 
respect to the coin Hilbert spaces. The partial tracing corresponds to coin 
tossing in an ordinary random walk, and results into a density matrix for 
the quantum system of the walker, which further may provide statistics of 
various quantum observables of the walk. 
Explicitly the unitarization of ARW reads 

^V^Qw) ^ ^ Pmi---PmkS mi ...S rnk Q w ( K S rni ...S mk ^ 

m 1 ,...,m k =± 

= Trf k {V® k Q® k ®Q w V® k ^), (62) 
and can be identified with a QRW. 
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Let us remark at this point that an equivalent decomposition of the 
unitary V® k would be 

k 

v m = Y[W t , where W t = (P+U)i (g)(S+)fc+i + (P-U)i (g)(5_) fe+1 , 

i=l 

(63) 

and the subindex denotes the position of the embedding of the respective 
operator into the A:- fold tensor product. In fact each of these operators 
Wi e End(Hf k H w ), provide a new decomposition of the CPTP map 
i.e £y — £w 1 w 2 ---w k , which is equivalent to a nonstationary QRW with 
k different unitary evolution operators empoyed in order to construct the 
k — th step. As a matter of fact this new decomposition helps to account 
for the type of quantum entanglement involved between coin and walker 
systems. Let us take the simplest k = 2 case, where V® 2 = W\W 2l with 

W x = (P+(g)l c (g)5 + + P_(g)l c (g)S_)^(g)l c (g)l 10) 

W 2 = (l c (g)P+(g)S + + l c (g)P_(g)S_) l c (g)U<g)l w - (64) 

The action of these operators on the product density matrices p c p c Pw 
is akin to theaction of some control-control-S'± type of non-local operator, 
which uses the two coin states as control spaces and the walker state as 
the target space, preceded by the local unitary operator U which acts on 
the control spaces and creates appropriate superposition of coin states. 
These actions generate quantum entanglement and can be described by the 
quantum circuit of Fig.l.c below. For purpose of comparison in Fig.l.a, we 
have included the corresponding circuit that generates the four entangled 
bipartite Bell states upon action of the composite operator UcnH <S> 1 on 
the four orthogonal product qubit states, and on Fig. lb the circuit that 
corresponds to the unitary V = (P+ S+ + P_ S-)U(g) 1, of the e v 
map. 

Though the topic of the entanglement in the QRWs will not be inves- 
tigated further here, it should be obvious from the above analysis that the 
CPTP map that implements some discrete time-step of a QRW, also gen- 
erates quantum entanglement that can be studied by appropriate circuits 
and evaluated by effective measures, as usually is done in other cases of 
coupled quantum systems. Although quantum correlations have been gen- 
erally accepted to be the common cause of all novel effects in the QRW 
performance, the exact evaluation of the entanglement resources needed in 
the course of a QRW is still an open problem. 
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Figure 1. Fig.l (a) Circuit generating entangled Bell states from product states using 
local Hadamard and non-local control-not gate; (b) Circuit generating the unitary evo- 
lution of the 1-step map ey ; (c) Circuit generating entangled coin-walker states from 
corresponding factorized ones using local unitary (e.g Hadamard gate in the case of the 
homonymous QRW), and non-local control-control- S± gates, that form the map of the 
2-stcp QRW e v = £WiW 2 - 



To establish further the connection among ARQs and QRWs we give 
four particular examples: 

i) the Euclidean QRW 13 [iso(2)-QRW]: with the distance operator 
L \m) = m \m) , with its eigenspace =span{|m) ,m £ Z}, and its dual 
phase operator $>\(p) — (f\(f) , with its eigenspace —span{\cp) , ip £ 
[0, 27r); o 2 }, related by a Fourier transform with H^. Two ARWs and its as- 
sociated QRWs (modulo local unitary operators in coin spaces, as explained 
above), can be constructed: the distance random walk on Z, with CPTP 
map constructed with Kraus generators been the step operators in the dis- 
tance operator eigenstates i.e S± \m) = E± \m) — e ±l * |m) = |m ± 1) , and 
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the phase random walk on the circle S, with CPTP map constructed with 
Kraus generator been the step operators in the phase operator eigenstates 
i.e S±\ip)=e tiL \ip±l); 

ii) the Canonical Algebra QRW [ hw-QRW]: with the position operator 
Q\q) = q \q) , with its eigenspace H® =span{|q) ,g e R; dq} 7 and its dual 
momentum operator P\p) — p \p) , with its eigenspace H p =span{ \p) , p S 
R; dp}, related by a Fourier transform with H®. Two ARWs and its asso- 
ciated QRWs (modulo local unitary operators in coin spaces, as explained 
above), can be constructed: the position random walk on R, with CPTP 
map constructed with Kraus generators been the step operators in the 
position operator eigenstates i.e S±\q) = e ±lP \q) = \q±l), and the 
momentum random walk on R, with CPTP map constructed with Kraus 
generator been the step operators in the momentum operator eigenstates 
i.e S ± \p) = e ±i Q\p) = \p±l) ] 

Hi) the M— dimensional Discrete Heisenberg Group QRW [Iim-QRW]: 
with the action operator N \ n) = n \ n) , with its eigenspace = span{|n > 
, n e Zm}, and its dual angle operator 8 |i? n ) = d n \d n ) , with its eigenspace 
H® =span{|^„) , e ^Zm}, related by a finite Fourier transform with 
H%. Two ARWs and its associated QRWs (modulo local unitary operators 
in coin spaces, as explained above), can be constructed: the action random 
walk on Zm , with CPTP map constructed with Kraus generators been the 
step operators in the action operator eigenstates i.e S± \n) = e ±ie \n) = 
/i* 1 \n) = \n ± 1) , and the angle random walk on j-Z^, with CPTP map 
constructed with Kraus generator been the step operators in the momentum 
operator eigenstates i.e S± \i} n ) = g ±x |i?„) = e ±i/N |i? n ) = |i? n ±i) ; 

iv) the Coherent State QRW 15 [ CS-QRW]: with the annihilation op- 
erator a\a) — a \a) , with its eigenspace =span{|a) , a e C; ^}.Two 
ARWs and its associated QRWs (modulo local unitary operators in coin 
spaces, as explained above), can be constructed: the annihilation ran- 
dom walk on C, with CPTP map constructed with Kraus generators been 
the step operators in the annihilation operator eigenstates i.e S±\a) = 
e^* 1 ) | a) = D±\ \a) = \a ± 1) . The step operators here identified as spe- 
cial case of the canonical coherent state displacement operator D±p = 
e ±/3at=F/3a = e J ( ±l3 \ have based the indicated step property on the following 
operator identity D a Dp = e~ lax/3 D a+ p, applied for the case of co-linear 
a, (3 vector on complex plane. 

To give an explicit identification of the ARW based on the hw alge- 
bra constructed in 15 , as a quantum random walk, and in particular as a 
CS-QRW, we make the following choices: the transition probabilities are 
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p+ = p,p- = 1 — p, the coin state is p c = |0 >< 0|, the U operator is 
U(p, 1 — p), and the step operators are CS displacement operators with 
steps ±/3 G C. Then the total 1-step evolution operator in the coin+ walker 

system is V = £ m _ ± P m U ® S m = ( vT ^f +/3 ) , and the 

reduced walker evolves in 1-step by the CPTP £v{Qw) = pD+pQ w D' + p + 
(l-p)D_^ft DL / 3=Tr c (^ c ®ft„Vt). 

For n steps the evolution of the walker has been chosen in 15 to be 
Sy(g w ). It is important to notice that this is a choice based on the ARW 
construction methodology, and that our present treatment of the same walk 
as a QRW, sees the e v {q w ) type of evolution to result from a partial tracing 
of the coin system at every step. Our previous discussion of other types of 
tracing schemes motivates the study of CS-QRWs with delayed tracing, in 
order to investigate phenomena such as enhanced or anomalous diffusion in 
ARWs. This problem will be taken up elsewhere. 

5. Discussion 

We have outlined a mathematical framework where the conception of ran- 
dom walk and its associated statistical notions, and equations of motion, 
can both be studied in an algebraic and quantum mechanical manner. 
ARWs and QRWs appear to be two aspects of the same mathematical 
device, so their interconnection serves to conceptually clarify the common 
ground between them and to enrich the heuristics of formulating new prob- 
lems and methodically searching for their solutions. 

Quantum random walks arc important both as quantum algorithms to 
experimentally be realized and as modules in a general quantum computing 
algorithm-devise that could outperform some classical rival. The connection 
ARW-QRW could serve to generalize, unify and compare such algorithms. 

Also Quantum Information Processing concepts and tools, could be de- 
veloped for ARW-QRWs. The step taken here is only a preliminary one 
towards developing such a theory. 

Finally, ARW-QRWs come with lots of free choices for its constituting 
parameters. To mention only one expected application in the field of Open 
Quantum Systems, we should emphasize the importance of choosing the 
functional in e.g the hw ARW. Various choices of functionals in terms of 
types of coherent state vectors, combined together with various choices of 
ordering the operator basis in the enveloping algebra i/(hw), i.e normal, 
antinormal, symmetric etc, could serve as a guiding rule for constructing 
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quantum master equations for open boson systems interacting with various 
types of quantum mechanical baths. 
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